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1 Introduction 

Given a multidimensional Brownian motion W on a probability space, consider the system of 
forward and backward stochastic differential equations 

(Xt =x + f bs{Xs, Ys) ds + f asdWs 

\Yt =h{XT) + If gs{Xs,Ys,Z,)ds- yf ZsdWs, t e [0,T] 

where x is a known initial value, T > 0 and b, a, g and h are given functions. In this paper, 
we give conditions under which the system admits a unique solution in the case where the value 
process Y is multidimensional and the generator g can grow arbitrarily fast in the control process 
Z. 

The focus in this paper is on Markovian systems, in which the functions b, a, g and h are deter¬ 
ministic. We consider generators that are Lipschitz continuous in X and Y and locally Lipschitz 
continuous in Z. In this setting, and for one-dimensional value processes, the decoupled system 
(with b depending only on X) has been solved by Cheridito and Nam [5] based on Malliavin 
calculus arguments. In fact, using that for Lipschitz continuous generators the trace of the Malli¬ 
avin derivative of the value process V is a version of the control process, they showed that the 
control process can be uniformly bounded, hence enabling solvability for locally Lipschitz gen¬ 
erators by a truncation argument. We make ample use of their method to deal with the coupled 
system (1.1). In this case, we propose a Picard iteration scheme which yields a sequence that can 
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be proved to be a Cauchy sequence in an appropriate Banach space under uniform boundedness 
of the control processes derived using Malliavin calculus arguments provided that the time hori¬ 
zon is small enough. Moreover using the PDE representation of Markovian Lipschitz FBSDEs 
as developed for instance in Delarue [7] and a pasting precedure, we construct a unique global 
solution for generators with growth specified on the diagonal and under additional assumptions, 
mainly non-degeneracy of the volatility a, see Theorem 2.2. 

Systems such as (1.1) naturally appear in numerous areas of applied mathematics including 
stochastic control and mathematical finance, see Yong and Zhou [28], El Karoui el al. [11] and 
Horsl el al. [16]. As shown for insfance in Ma el al. [23] and Cheridilo and Nam [5], in the 
Markovian case, EBSDEs can be linked to parabolic PDEs. More recently Fromm et al. [14] 
proved that FBSDEs can be used in the study of the Skorokhod embedding problem. 

BSDEs and FBSDEs with Eipschitz continuous generators are well understood, we refer to 
El Karoui et al. [11] and Delarue [7]. If Y is one-dimensional and g is allowed to have quadratic 
growth in the control process Z, BSDEs’ solutions have been obtained by Kobylanski [20], 
Barrieu and El Karoui [2] and Briand and Hu [3, 4] under various assumptions on the terminal 
condition ^ = h{XT)- We further refer to Delbaen et al. [8], Drapeau et al. [10], Cheridito 
and Nam [5] and Heyne et al. [15] for results on one-dimensional BSDEs and FBSDEs with 
superquadratic growth. Mainly due to the absence of comparison principle, general solvability 
of multidimensional BSDEs with quadratic growth is less well understood. Some important 
progress have been achieved recently for BSDEs with small terminal conditions, see Tevzadze 
[26] and for more recent development, see Hu and Tang [17], Euo and Tangpi [22], Jamneshan 
et al. [19], Cheridito and Nam [6], Frei [12] and Xing and Zitkovic [27]. 

To the best of our knowledge, the only works studying well-posedness of coupled FBSDEs with 
quadratic growth are the article of Antonelli and Hamadene [1] and the preprints of Euo and 
Tangpi [22] and Fromm and Imkeller [13]. 

In [1] the authors consider a one-dimensional equation with one dimensional Brownian motion 
and impose monotonicity conditions on the coefficient so that comparison principles for SDEs 
and BSDEs can be applied. A (non-necessarily unique) solution is then obtained by monotone 
convergence of an iterative scheme. In [13], a fully coupled Markovian FBSDE is considered 
with multidimensional forward and value processes and locally Eipschitz generator in {Y, Z) 
and a existence of a unique local solution is obtained using the technique of decoupling fields 
and an exfension fo global solufions is proposed. Allhough the (non-Markovian) system studied 
in Euo and Tangpi [22] is the same as the one considered in the present paper, the techniques 
here are essentially different. Furthermore, the main results we present here can be extended to 
the non-Markovian setting and to random diffusion coefficient (when a depends on X and Y) 
under stronger assumptions involving the Malliavin derivatives of g and h. We refer to the Ph.D. 
thesis of Euo [21] for details. 

In the next section, we present our probabilistic setting and the principal results of the paper. In 
Section 3 we prove local solvability of multidimensional BSDE with superquadratic growth and 
give conditions guaranteeing global solvability. Section 4 is dedicated to the proofs of the main 
results. 
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2 Main results 


Let be a filtered probability space, where {J^t)te[o,T] is the augmented 

filtration generated by a d-dimensional Brownian motion W, and T = Tt for a finite time 
horizon T E (0, oo). The product x [0, T] is endowed with the predictable (i-algebra. Subsets 
ofM*^andR^''^A;EN, are always endowed with the Borel (i-algebra induced by the Euclidean 
norm | • |. The interval [0,T] is equipped with the Lebesgue measure. Unless otherwise stated, 
all equalities and inequalities between random variables and processes will be understood in 
the P-almost sure and P (g) df-almost sure sense, respectively. For p E [l,oo] and A: E N, 
denote by the space of all predictable continuous processes X with values in such 

that := P[(supjg[o,T] l^tl)^] < by the space of all predictable 

processes Z with values in such that := E[{J^ < oo. 

Let l,m E N be fixed. The solution of (1.1) with values in R"^ x R^ x R*^'^ can be obtained 
under the following conditions: 

(Al) 6 : [0, r] X R^ —)• R"* is a continuous function and there exist ki,k 2 , Xi > 0 such that 
\bt{x,y) - bt{x',y')\ <ki\x- x'\ + A;2 Id - y'\ and \bt{x,y)\ < Ai(l + |x| + |y|) 
for all x, x' E R”^ and y, y' E R^ 

(A2) a : [0,T] ^ measurable function and there is A 2 > 0 such that \at\ < A 2 

for all t E [0, T]. 

(A3) h : R™ —)• RMs a continuous function and there exists /cs > 0 such that 

\h{x) — h{x')\ <k^\x — x'\ 

for all X, x' E R”^. 

(A4) g : [0, T] x R™ x R^ x R^^®* —)• RMs a continuous function satisfying y.(0,0,0) E 
L^([0, T]), and there exist /cs, /C 4 > 0 and a nondecreasing function p : R+ —)• R+ such 
that 


\gt{x,y,z) - gt{x\y,z)\ < A :3 |a: - x'| ( 2 . 1 ) 

for all X, x' E R”^, y E R^ and z E R^^'^ such that \z\ < M := AX 2 k^'/dl and 
|yt(x,y,z) - gt{x,y',z')\ < A: 4 |y-d'| + p{\z\ V |/|) 
for all X E R™, y, y' E R^ and z, z' E 
(A5) There exists a constant K > 0 such that 

\gtix,y,z) - gt{x',y,z) - gt{x,y',z') + yi(x', y', z')| <K\x- x'| (|y - y'| + \z - z'\) 
for all t E [0, T], x, x' E R”", y, y' E R' and z, z' E 
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Our first main result ensures local existence and uniqueness for the coupled FBSDE (1.1) under 
the previous assumptions. The proof is given in Section 4. 

Theorem 2.1. Assume that (A1)-(A5) hold. Then there exists a constant Ck^x^i^d > 0 depend¬ 
ing only on ki, k 2 , k^, k^, k^, A 2 , I and d, such that the FBSDE (1.1) has a unique solution 
{X,Y,Z) G 52(M”") X 52(M') X with\Zt\ < M, whenever T < CkxiA 

Local existence results as Theorem 2.1 above have been obtained in [13, Theorem 3] and [22, 
Theorem 2.1] in essentially different settings and with different methods. A natural question is to 
find conditions under which the above result can be extended to obtain global solvability. Fromm 
and T mk eller [13] propose a concatenation procedure allowing to prove existence of solutions 
to a fully coupled FBSDE on a "maximal interval". In the present setting, under additional 
assumptions, a pasting method based on PDEs allows to get global existence and uniqueness for 
the FBSDE (1.1). The proof is also given in Section 4. 

(A4’) g : [0,T] x x x ^ is a continuous function satisfying g.(0,0,0) G 
L^([0, T]), and there exist ks, k 4 > 0 and a nondecreasing function p : M+ —)• M+ such 
that 

x'l + fc4 ly - g'l + p[\z\ V \z'\) |/ - (z')*| 

for alH = 1,..., X, x' G M™, y, y' G and z, z' G 

Theorem 2.2. Assume that (Al), (A2), (A3), (A4’), and (A5) hold and there exist A 3 > 0; 
A 4 , As > 0 such that^ 

\h{x,y)\ <Ai(l + |y|) 

{x,(Jtatx) > Aalxp 2) 

\gt(x,y,z)\ < A 4(1 + |y|+/ 5 (| 2 ;|) | 2 ;|) 

|/i(x)| < As 

for ail t G [0, T], x, x' G M™, y, y' G and z, z' G then the FBSDE (1.1) has a unique 

global solution (X,Y,Z) G 5^(]R™') x 5^(M^) x such that \Zt\ < M for some 

constant M > 0. 


Remark 2.3. The following counterexample shows that the additional conditions in Theorem 2.2 
cannot be dropped without violating global solvability. Consider the FBSDE 

(Xt = f/, Yu du 

kXu du - Zu dWu. 

This equation can be rewritten as 


Y = 


IS 1 

J j kYu duds — J 

t 0 t 


Zu dWn 


'(j( is the transpose matrix of at. 


(2.3) 
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It has been shown in [9, Example 3.2] that if T'/k < | then the BSDE with time-delayed 
generator (2.3) has a unique solution whereas if Ty/k = (2.3) may not have any solutions 

and if it does have one, there are infinitely many. ♦ 

Theorem 2.2 extends to fully coupled generators provided that a more special structure is as¬ 
sumed. The proof of the next proposition is given in Section 4.3. 

Proposition 2.4. If there exist an invertible matrix T E and A 3 > 0; A 4 , A 5 > 0 such that 
(Al), (A2), (A3), (A4’ ), (A5) and (2.2) hold for h := Th, gt{x, U, z) := T^^y, T^^z) 

and bt(x,y) := bt(x,r~^y), then the FBSDE ( 1 . 1 ) has a unique global solution {X,Y, Z) E 
5^(M™') X 5^(M^) X 5°°(]R^^‘^) such that \Zt\ < M for some constant M > 0. 

The point of the above proposition is that global solvability can be obtained for a generator 
that does not satisfy the "diagonally superquadratic" growth condition (A4’) provided that the 
transformed generator g does. Moreover, notice that Theorem 2.2 and Proposition 2.4 guarantee 
existence of a decoupling field, see [13]. In particular, the boundedness of Z guarantees a 
uniform Eipschitz property of the decoupling field. 

Theorem 2.1 relies on an existence result for multidimensional BSDEs presented in Nam [24] 
and revisited in the next section. 


3 Multidimensional BSDEs with bounded Malliavin derivatives 


Eet us introduce the spaces of Malliavin differentiable random variables and stochastic processes 
'D^’^(M^) and ^(R^). Eor a thorough treatment of the theory of Malliavin calculus we refer to 
Nualart [25]. Eet M. be the class of smooth random variables ^ = (^^,..., ^^) of the form 


T T 

e = I hTdWs) 

0 0 

where yp® is in the space C'“(R”;R) of infinitely continuously differentiable functions whose 
partial derivatives have polynomial growth, ..., /i*®® E E^([0, T]; R®^) and n > 1. Eor every 
^ in A4 let the operator D = (D ^,..., D'^) ■. M. ^ Lf{Vl x [0, T]; R®^) be given by 



0 < f < T, l<i<l, 


and the norm ||^||^ 2 ■= + Jo As shown in Nualart [25], the operator D 

extends to the closure P^’^(R^) oftheset A4 with respect to the norm H-H^ 2 - A random variable ^ 
is Malliavin differentiable if ^ E P^’^(R^) and we denote by Dti its Malliavin derivative. Denote 
by £a^(R^) the space of processes Y E 4f^(R^) such that T) E P^’^(R^) for all t E [0,T], the 
process DYt admits a square integrable progressively measurable version and 



T T 


+ E 


\DrYff dr dt 


< 00. 


0 0 
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We next consider a system of superquadratic BSDEs of the form 

T T 

Yt = C + J gu(Yu, Zu)du - J ZudWu (3.1) 

t t 

satisfying the following conditions: 

(Bl) g : n X [0,T] x x —)• RMs a measurable function and there exist a constant 

B G M+ and a nondecreasing function p : R+ —)• R+ such that 

\gt{y, z) - gt{y', z')\ <B\y-y'\ +p[\z\y \z'\) \z - z'\ 
for all t G [0, T], y, y' G R* and z, z' G 

(B2) ^ G Il^’^(R^) and there exist constants Aij > 0 such that \Di^^\ < Aij for all i = 
1,... ,1, j = 1,... ,d and t G [0,T], 

(B3) 5 .(0,0) G 2(^(R^) and there exist Borel-measurable functions qij : [0,T] —^ R+ 
satisfying qfj{t)dt < 00 and for every pair (y,z) G R^ x R^^'^ with 


z\ < Q ■■= 



it holds 

• g.iy,z) G £a^(R') with \Digl{y,z)\ < qij{t) for alH = 1,j = 1,d 
and u G [0, T], 

• for almost all u £ [0, T] one has 


Dugt{y,z) - Dugt{y',z')\ < Ku{t) (|y - y'\ + \z - z'\) 


for all f G [0, T], y, y' G R* and z,z' G R^^'^ for some R_|_-valued adapted process 
{Ku{t))t(^[Q,T] satisfying ||i^«||^ 4 (]R) du < 00 . 

The following is an extension of Cheridito and Nam [5, Theorem 2.2] to the multidimensional 
case. It was proved in Nam [24] under slightly different assumptions. We give the proof for the 
sake of completeness. 

Theorem 3.1. Assume that (B1)-(B3) hold and T < Then the BSDE (3.1) admits 

a unique solution in 5^(R^) x 5°°(R^^'^) and \ Zt\ <Q. 


Consider the following stronger versions of the conditions (B1) and (B3): 

(B1’) y is continuously differentiable in (y, z) and there exist constants B G R+ and p G R+ 
such that \dygt{y,z)\ < B and \dzgt{y,z)\ < p for all t G [0,r], y,y' G R' and 
z,z' G 
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(B3’) The condition (B3) holds for all (y, z) G x 


Lemma 3.2. ^fBl’j, fB2j and fB3’j hold, then the BSDE (3.1) admits a unique solution 
(Y,Z) G X and 



(3.2) 


Proof. By Cheridito and Nam [5, Lemma 2.5], the condition (B2) implies E [|^|p] < +oo for all 
p G [1, oo). It follows from El Karoui et al. [11, Theorem 5.1 and Proposition 5.3] that the BSDE 
(3.1) has a unique solution (V, Z) G 5^(M^) x 2(^(]R^^'^), which is Malliavin differentiable. 
Moreover for every i = 1,... ,l and j = 1,... ,d, the process {DrY^, DrZ^*)ie[o,T] has a 
version )te[o,r] which satisfies 

= 0, = 0, for 0<t<r<T, 


and is the unique solution in 5^(M*) x of the BSDE 

T T 

ut = DiC + j dydsiYs, Zs)W/ + d.gsiYs, Z,)L/> + D^y.(y„ Z,)ds - J 


- / Vf^dW.,. 


Applying Ito’s formula to yields 


C //>|2 = | Di ^|2 _ / 2 W / Vf^dWs 


+ I 2W/dygs{Ys, Zs)Uf^ + 2W/d,gs{Ys, Z^W^ + 2Uf^Dlgs{Ys, Zs) - 


<\Dlif - / 2W/Vf'^dWs+ / 2B\W/f+ 2p\Uf'^\\Vf^\+2 


< 


V 

T T I 

mi?-j 2u>/vmdw, +1 {2B+pmi)iw/? + J2qHs)d, 


EiSwit/ri-iv/’n"* 


t t 

Using condition (B3) and taking conditional expectation in the above inequality yields 

T 


\ui’r < + [{‘^B + p^ + 1) + ^qUs)ds 


i=l 


2=1 


Tt 


(3.3) 


By El Karoui et al. [11, Proposition 5.3] the process Z is a version of the trace (U/)ig[o,r] of the 
Malliavin derivative of Y. Hence (3.2) follows from (3.3) by applying GronwalTs inequality. □ 
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Proof (Theorem 3.1 ). Define the Lipschitz continuous function g by 


9t{y,z) 


9tiy,z) if \z\<Q, 

9 t{y,Qz/\z\) if \z\>Q. 


(3.4) 


By Cheridito and Nam [5, Lemma 2.5] and El Karoui et al. [11, Theorem 5.1] the BSDE corre¬ 
sponding to (g, f) has a unique solution (Y, Z) G x Eor x = (y, z) G ]^i+i-xd. 

let /3 G be the mollifier 

if W<1. 

[O otherwise, 


where the constant A G M+ is chosen such that /jjj+ixd I3{x)dx = 1. Set /3'^{x) := v}~^^^‘^(5(nx), 
n G N \ {0}, and define 

g^{uj,x):= J gt{uj,x')/3'^(x - x')dx' 

Ri + ixd 


SO that for each n > 0 the function g"^ satisfies (BE) and (B3’) with the constant p replaced 
by p{Q). By Eemma 3.2 the BSDE corresponding to ( 5 ”, has a unique solution {Y^, Z^) in 
X which satisfies 


i=i i=i 

e(2B+p^(Q)+l)T^ 

Since T < 

for all j = 1 ,..., d. 

This shows \Zf\ < Q. Since p" converges uniformly in (t, ui, y, z) to g, using the procedure of 
the proof of Cheridito and Nam [5, Theorem 2.2], it follows that (Y'^, Z'^) converges to (Y, Z) 
in X so that \Zt\ < Q. Since g{y,z) = g{y,z) for all {y,z) G x 

with \z\ < Q, it follows that (Y, Z) is the unique solution of the BSDE corresponding to (^, g) 
in 5 ^(M') X □ 




j g( 2 B+p 2 (Q)+l)(T-t) 


Corollary 3.3. Suppose fBlj-fBS) hold, T < 2 b+2h2+^ ^ xS°°{R^^^) is 

the solution of the BSDE (3.1). Then Yt G D^’‘^(R^)for all t G [0, T] and for every j = 1,... ,d, 
one has 


\DiYtf<2(j2^h + i2 I 


2=1 


2=1 


for all r G [0, t]. 


(3.5) 





Proof. Since \Z\ < Q is bounded, {Y, Z) solves the BSDE with terminal condition ^ and gen¬ 
erator g defined by (3.4). If g satisfies (Bl’) and (B3’), then the result follows from Lemma 
3.2. Otherwise consider the sequence of smooth functions g'^ converging to g as defined in 
fhe proof of Theorem 3.1. Let (Y^,Z^} G 5^(]R^) x be the solutions to the BS- 

DEs corresponding to ( 5 ”,^), which converge to {Y,Z) in x By Lemma 

3.2 {Yf,Zf) G X for each t G [0, T] and the arguments in the proof of 

Theorem 3.1 imply 

qfj{s)d^ j = I,..., d, r, f G [0, T]. 

Hence, sup^^pj Elf^ iD^Yfl'^ dr] < 00 for each t G [0, T]. Since (Yf) converges to 1) in 
it follows from Nualart [25, Lemma 1.2.3] that Y) G and (DYf) converges to DYt in 

the weak topology of Thus, DrYt satisfies (3.5). □ 



As a concequence to Theorem 3.1, we give a condition for global solvability of fully coupled 
systems of BSDEs. Lor the remainder of this section we put 


log(2) 

2B + p‘^i2^Q) + V 


n G N. 


Proposition 3.4. Assume that (B1)-(B2) hold, that there exists G N such that Yln=o > 
T, and (B3) holds with Q replaced by 2^Q. Then the BSDE (3.1) has a unique solution in 
X and\Zt\ < 2^Q. 


Proof. If T < Aq then the result follows from Theorem 3.1. Otherwise, if T > Aq it follows 
by the same arguments as in the proof of Theorem 3.1 that the BSDE (3.1) has a unique solution 
{Y^, Z^) in 5^(M^) x S°°{R^^^) on the interval [T — Aq, T]. Moreover, Z^ satisfies \Z^ \ < Q 
and by Corollary 3.3 one has G P^’^(M^) and for every r <T — Aq, 


ID^Y^-Ao 


2 



Since g satisfies (B3) for all {y,z) G x such that \z\ < cQ, again by Theorem 3.1 
the BSDE (3.1) with terminal condition Yrp_^^ has a unique solution (Y^,Z^) in 5^(]R^) x 
^oo(]^Zxd) QJJ [(j. - Ao - Ai) V 0, T - Ao] , and 


l^^>^(T-Ao-AOvol' + E / (2' + 


i=l 


2=1 


\Zl\ < 2Q, t G [(T — Aq — Ai) V 0, T — Aq] 


for all j = 1 ,..., d 


Repeating the previous arguments, for A > 2 the BSDE (3.1) has a unique solution {Y^, Z^) 
in 5^(]R') X S°°iR^^'^) on [(T - A„) V 0, (T - A„) V 0] with terminal condition 
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^(r-E^L-o^A.)vo-Moreover, 


l T N 




n=0 


n=0 


+ forallj = l,...,d 

i=l i=l Q k=l 

N N-l 

|Zf| <2^Q, t G [(r-^A„) V0,(r- ^ A„) VO . 

T) 

Hence, the pair (E, Z) given by 

Y := E°l[r-Ai,T] + X] ^"^[(t-ELo A*)vo,(r-EILV Ai)vo] 

n=l 
N 

Z := E°l[r_Ai,T] + Y1 ^’"^[(r-E?=o Ai)vo,(T-Eiroi Ai)vo] 


N 


n=l 


solves (3.1) and its uniqueness follows from Theorem 3.1. 


□ 


Remark 3.5. The condition > T for some N £ N does not guarantee global solv¬ 

ability of multidimentional BSDEs with superquadratic growth. In fact, if p{x) > C(1 + y/x) 
for all X > 0, then X]n>o However, it does guarantee global solvability for BS¬ 

DEs whose generator grows slightly faster than the linear function. Eor instance, if p{x) < 
(7(1 -h Y^log(l -h x)) one has 

log(2) 


E 


n =0 
oo 

E 


log(2) 


2B + p2(2^Q) + 1 - ^ 2H + 2(72(1 + log(2A^(l + Q))) + 1 

log(2) 


2B + 2(72(1 -h log(l -h Q) -h nlog(2)) -f 1 


= oo. 


4 Coupled FBSDE with superquadratic growth 

4.1 Proof of Theorem 2.1 

Step i: We first assume that h, b and g are continuously differentiable in all variables. Eet us 
define 

^1 M A A A log 2 

kl ki k2M 2A:4 + p2(M) + l 

with M := 4^5A 2 \/^. We will show that for T < the sequence {X"^, Y^, Z^) given 

by = 0, E° = 0, ZO = 0 and 

f Af+i = X + £ 6(A-+i, E-) du + £ CT, dWu 
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is well defined and that \Z^\ < M for all n G N and t € [0, T], The process is well defined, 
Xl belongs to for every t and the process (T*r-^t)te[o,r] satisfies the linear equation 

DrXl = 0, 0<t<r<T, 

au , 0 < r < t < r, 

with Dr{f^ au dWu) = <7l[r,^]^ see Nualart [25, Lemma 2.2.1 and Theorem 2.2.1]. Hence, since 
b is Lipschitz continuous, we have 


DrX} = / {d^bDrXl + dybDrY^) du + Dr 


iDrXU < / kiDrXldu + ar 


and 


DrXl \ < X2e^^\ 


where the second estimate comes from Gronwall’s inequality. We will now show that since T < 
XI d’ h{X)p) and g{X^, •, •) satisfy (B1)-(B3). In fact, since h is continuously differentiable 
and X}p G it follows from the chain rule, see for instance Nualart [25, Proposition 

1.2.4], that h{X^) G and \Dl{h{X^))\ = \d^h{X^)Dix^\ < for all r G 

[0, T], j = 1,... , d. Using T < we deduce that h{X}p) satisfies (B2) with Aij := 2X2h^- 
Similarly, by (A4) and using that the function x ^ g{x, y, z) is continuously differentiable, 
it follows that gXX^,y,z) G and \D^{g]{X^^y, z))\ < A 2 A: 3 e^^L j = 1,... ,d for 

all {y-iZ) G X such that \z\ < M and, due to (A5), applying the same argument to 
gt{x, y, y\ z, z') := gtix, y, z) - gt{x, y\ z') yields 

\Dlgt{X^,y,z) - D^gtiX},y', z')\ < KX2e^^\ 


Using T < || A we deduce that g, {X^,y, z) satisfies (B3) with qij = 2X2kz and Ku{t) := 

2KX2- Moreover due to (A4), the function (t, y, z) i-A gt{Xl,y, z) satisfies (B1). 

Therefore, by T < 2 k 4 +p^M)+i ’ Theorem 3.1 ensures that (U^, Z^) exists. Consider the func¬ 
tion g defined by 

- . . i 9 t{x,y,z) if|z|<M 

qtix, y,z) = < 

\gt{x,y,zM/\z\) if \z\>M. 

Since (U^, Z^) also solves the BSDE with terminal condition h{Xi^) and a Lipschitz generator 
g{X^, ■, •), it follows from Lemma 3.2 and its proof that (1^^, Z^) G x for 

all t G [0,T] and DtY^ is bounded and it holds Zil = DtY^. In addition, we have < 

4A2 and |L>r-A"/| < M. 

Now let n G N, assume that {X^,Yt^, ) G x x Zf = 

and |AA”| < 4 A 2 , \DrYl^\ < M for all r, t G [0,T]. The process is well defined, for 
each t; X^~^^ belongs to and it holds 

DrX^+^ = 0, 0 < t < r < T, 
t 

DrX^^^ = (^r + I (A6AA:”+^ + dybDrY^) du, 0 < r < t < T. 
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Since dxb, dyb and a are bounded by ki, k 2 and A 2 respectively, it follows from Gronwall’s 
inequality that 

\DrX]}+^\ < + A :2 J \DrY^\ du ) . 

\DrX^+^\ < (A2 + k2TM) < 00 


Hence, 


(4.1) 


so that since T < < 4 A 2 . As above, and p. y, z) are 

Malliavin differentiable and satisfy (B1)-(B3) with Aij := 2X2k^, qij = 2X2kz and Ku{t) := 
2KX2- It then follows again from Theorem 3.1 that exists and < M 

is bounded. Since also solves the BSDE with terminal condition h{X^'^^) and 

a Lipschitz generator •, •), Lemma 3.2 and its proof guarantee that (1^”“'"^, E 

X pA 2 (]^zx<i) fQj. ^ g [Q DtY^+^ is bounded and it holds Xf+i = 

with < M. 

Step 2: Now we show that there is a positive constant 2| ^ ^ ^ such that if T < 2\xi , 1 , then 
is a Cauchy sequence in S^{W^) x x Using (Alj we can 

estimate the norm of the difference Xf — Xf as 


|X”+^ - Xf < 2 




Thus 


sup |Xf+^-Xf|2<2 
0<t<T 




Taking expectation on both sides and using Cauchy-Schwarz’ inequality, we have 





- rp 


- J, 

E 

sup ixr+i-xr|2 

< 2TklE 

[ |X"+1 -X^pds 

+ 2TklE 

[ |y”-y”-ipds 


0<t<T 

J 

.0 

J 

.0 


< 2T‘^klE 

sup 

+ 2T‘^klE 

sup ly.’^-y^-'i^ 


0<t<T 


0<t<T 


Choosing T to be small enough so that 2T‘^kf < it follows 


E 

sup |xr+i-xr|2 

< AT'^klE 

sup |y--y--Y 


0<t<T 


0<t<T 
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On the other hand, applying Ito’s formula to /3 > 0, we have 

T 

g/3t|yn+l _ yn|2 ^ _ /j(x")|2 _ 2^ e^\Y^+^ - - Z2)dWs 

t 

T T 

- j - Z'^fds - J - Yj^fds 

t t 

T 

+ 2 I - YD [gs{X2+\Yr+\ Z^^) - g,{X2, Y^, Z^] ds. 

t 

Hence, due to the condition (A3) and the boundedness of (Z*^), it holds 

T 

g/3t|yn+l _yn|2 ^ 


;”+^ - zD'^ds 


< e 


I e^^{z: 

t 

T 

\h{XD^) - h{XD\^ -2 J - y,’")(y"+^ - Z^dWs 

t 

T T 

- J - YD^ds + 2 j e^D{M) |y"+i - y,'^] \ZD^ - Z2\ ds 

t t 

T T 

+ 2 y |y”+i - y”| |X"+1 -x^\ds + 2 J |y”+^ -y"|^ ds. 

t t 

With some positive constants ai, a 2 , it follows from (A3) and Young’s inequality that 

T 


g/3t|yn+l _ y«|2 y 


- 2 


J - ZD'^ds < e^^kl\XD^ - X 

t 

T T 

D‘iYr^-Y-){zr^-zr)dw. + 012 I g/ 3 .|yn+l 


n\2 

t\ 


-xri^ds 


+ 


({p{M)r , kl 

V at 


+ ^ + 2/^4 - /3^ J e^"(y”+i - YP)^ds + aij - Z^\Ds. 


(4.3) 
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Letting /3 = —\- -^ + 2ks and taking expectation on both sides above, we have 


E 


■t 


2 +E 

j - Z2fds 

< e^^klE 


- j. 


■ T 

+ aiE 

f 

+ a2E 



J 


J 


Putting ai = I and a 2 = 1, the previous estimate yields 


E 

■ T 

f e/3*(^n+l _ 

< 2e^^kjE [|X”+^ - + 2E 

■ T 

[ e^"|X,^+i -X^pds 


J 

Lo J 


J 

Lo J 


Next, taking conditional expectation with respect to Ft in (4.3) gives 

T 


g/3t|yn+l + E 




Ft 


< e^^klE - XT^\^\Ft] 



■ T 



■ T 


+ aiE 

J e^^\Z^+^ - Z^fds 

.t 

Xt 

+ 02-E 

J e^"|X"+i - X.^l^ds 

.t 

Xt 


Thus, by Burkholder-Davis-Gundy’s inequality, with a positive constant ci and ai = ^, 02 = 1 
we have 


E 


sup -Y^\^ 


0<t<T 


< cief^^klE [|X”+i -X^|2] 


+ C1-E 


< 2cieP'^klE - X^\^] + 2ciE 

It now follows from (4.2) that 


■ T 


- rp 

f f,F\z^+i _ z^\^ds 

+ Cli? 

[ e^"|X,”+i -X^pds 

J 

Lo J 


J 

Lo J 


eF\x^+i - X^\^ds 





- J, 

E 

sup |y”+i-y"|2 

+ E 

[ (Z”+' - Z^fds 


0<t<T 


J 




.0 


< 8(ci + l)e^^(fe| + T)T^klE 


sup \Yr-Y,^-^\^ 
0<t<T 
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Taking T small enough so that 


8{ci + l)el^^{kl+T)T^kl<^, 

we obtain that {X^, Y^, Z^) is a Cauchy sequence in S'^{W^) x x Thus, it 

suffices to define 2| ^ ^ by fhe condifions 

12T^kl < i 

\s{ci + l)e^^{kl + T)T‘^kl < i. 

By confinuify of 6, g and h we have fhe exisfence of a solufion (X, X, Z) in x52(M') X 

-^ 2 (]^/xd) qj: fBSDE (1.1) and if follows from fhe boundedness of (Z”) fhat \Zt\ < M. The 
uniqueness in x x follows from fhe boundedness of Z and by re- 

peafing fhe above argumenfs on fhe difference of fwo solutions. 

Step 3\ If one of fhe funcfions b, g ox h is nof differentiable, we apply fhe fechnique of fhe 
proof of Theorem 3.1. Namely, we use an approximation by the smooth functions defined as 
follows: For n E N, lef /3^, and /3^ be nonnegafive funcfions wifh support on {x E : 
|a^| < {x E : |x| < and {x E ]R"*+N«xd . |^| ^ i| respectively, and satisfying 

(3i{r)dr = 1, I3n{r)dr = 1 and /g^+i+ixd I3n{r)dr = 1. We define fhe convolutions 

b7ix,y) ■■= j bt{x',y')/3l{x'- x,y'- y)dx'dy', h^{x) := j h{x')/3l{x'- x)dx', 

Rm+i Rm 

g'^iu, x, y, z) := j g{u, x',y', z')/3^{x' - x,y' - y, z' - z)dx'dy'dz'. 

+ Z X d 

If is easy fo check fhaf 6” satisfies (Al) wifh fhe consfanfs ki, k 2 and 2Ai and fhaf g"' and h"' 
safisfy (A4) - (A5) and (A3), respectively, wifh fhe same consfanfs. From Steps 1 and 2, fhere 
exisfs a posifive consfanf Ck^x,i,d independenf of n such fhaf if T < Ck,x,i,d, FBSDE (1.1) 
wifh paramefers {b'^,h"‘ , g"-) admits a unique solution {X^,Y'^,Z'^) E 5^(M”*) x 5^(M^) x 

5 oo(]^Zxd) 

\Zt\ < M. 

By the Fipschitz continuity conditions on b and h and the locally Fipschitz condition of g, the 
sequences (6”) and (/i”) converge uniformly to b and h on and M"*, respectively, and (g^) 
converges to g uniformly on x A for any compact subset A of Combining these 

uniform convergences with the boundedness of Z”, similar to above, we can show that there 
exists a constant Ck,x,i,d depending only on ki,k 2 ,ks,k 4 ,k^, X 2 ,l,d such that if T < Ck,x,i,d, 
(X^, y"', Z"’) is a Cauchy sequence in the Banach space 5^(R”^) x 5^(R^) x 4f^(R^^'^). 

In fact, for any m, n E N, using Cauchy-Schwarz’ inequality we have 

T 

\x--xn"<T I 
0 
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Thus, taking the supremum with respect to t and then expectation on both sides give 


\X^ - 


l52(R"*) 


< 3T 


1 

j y”) - y”)|" + KixT.rn - YT)i 


0 


T 

<3T [mx:,Y:) - b^{x:,Y:)f + i6:r(xr,>r) - 6„(x-,y-)|2) du 


+ 3klT‘^ IIX*" - X 


m\\2 


,-^+3k^T^ ||yn _ym||2 


(4.4) 


where the second inequality follows from (Al). On the other hand, applying Ito’s formula as in 
Step 2, one has 


|ym_yn| 2 ^ / 


1 

h 


n _ ym\2 

u I 


du 


< |/i’"(X^) - /i(XT)r + \h'^{X^) - h{XJp)Y + kl |X^ - XJp 

T 


- 2 


I {Y,^+^ -Yr){Z2+^ - Z2)dWs 


t 

T 


+ 1 \Y:-Y:r\{Kix:,Y:,z:)-g^{x:,Y:,z:)\ 


+ ^ ix- - x™i+ h\ y: - y^i 

+ p(M)|z:-y-|). (4.5) 

Taking expectation, due to Young’s inequality we have 


IZ*" - Z 


m\\2 


xd\ 


rm\\2 


< E[\h^iX^) - hixm + E[\h^iX^) - hiX^r] + (ki + -Tki) ||X" - X’^ 


1 of 

X ^ II, rm\\2 


_L " T IWYt vmu'^ I " / v"- r, ('yn\\‘i 

+ Ik ll52(R*) + ^ / \9u\^uXu ^^u)\ 


1 


+ \Tklp\M) ||y- - ^ \\z^ - 


On the other hand, taking conditional expectation in (4.5) and then the supremum with respect 
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to t and then expectation on both sides, we have due to Young’s inequality 


vrnu2 

^ ~ ^ ll52(]R*) 

< E[\h^{X^) - hiX^)\^] + E[\h^iX^) - hiX^)\^] + kl ||X" - 


+ -T I 
2 ' 


yn _ yTi 


.) + - / \g:{x:,Y: 


zr: 


-gu{x:,Y: 


zr: 


+ YT, - gu(x:r,Yr, ^r)r du 


1 


Tkl 


IX^ - X" 


1 


52 (Rm) + -Tk^p (M) \\Y 


-Y^‘ 


+ 


1 


2 o iio“(^iK""j ' 2 II- - lls^(R*) ' 2 

Combining (4.4) and (4.1) we observe that if T is small enough so that 


- Z" 


||2 

IIh2(k*x‘^) • 


13kfT^ < i 

I iT + 3klklT^ + iT^klkl + \Tklp‘^{M) < | 

then, the uniform convergence of ( 6 ”), g'^ and (/i”) to b, g and h ensure that {X^,Y^,Z^) 
is a Cauchy sequence. The verification that the limit (X, Y,Z) of the sequence (X"^, Y", Z”) 
solves the FBSDE (1.1) uses continuity of the functions b, h and g, and that \Zt\ < M is a conse¬ 
quence of the boundedness of {Z"^). Taking Ck,x,i,d ■= Ck,x,i,d/^Ck,x,i,d concludes the proof. □ 


Due to Theorem 2.1 above, our global existence result now follows from a pasting procedure. 


4.2 Proof of Theorem 2.2 


If T < Ck,x,i,d^ then the result follows from Theorem 2.1. 

Assume T > Ck^x,i,d and let hu : K —)• M be a continuously differentiable function whose 
derivative is bounded by 1 and such that h'j^{a) = 1 for all —M < a < M and 


fiMia) 


'(M + 1) 

< a 

-{M + 1) 


if a > M + 2 
if |a| < M 
if a < — (M + 2). 


An example of such a function is given by 


I ( \ - I “ ^)) ® ^ ^ + 2 ] 

^ ~ |(M2 + 2Ma + a(a + 4))/4 if [-(M + 2),-M], 

see Imkeller and Reis [18]. By the assumptions (A3) the function g : [0, T] x ME x ^ 
M defined by 

gt{x,y,z) := gt{x,y,hM{z)) 




(4.6) 


wifh hM{z) := {hM{z^^))ij is Lipschifz confinuous in all variables. Thus, if follows from 
Delarue [7, Theorem 2.6] fhaf fhe equation 


Xt=x + /q bu{Xu,Yu) du + /g cju dWu 

Yt = h{XT) + ~gu{Xu, Yu, Zu) du - dWu, t G [0, T] 
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admits a unique solution {X, Y, Z) € S'^{MX) x x Moreover, there exists 

a Lipschitz continuous function 6 : [0, T] x M*" bounded by a constant K such that 

Yt = 9{t, Xt) for all t G [0, T]. In fact, for every x, x' G M'^, t G [0, T] and i = 1,..., / we have 


e{t,x^)-e{t,xf) 

= h\X^) - h\X^) + 

= h\X?p) - h\x^) + 


„^(VX yx\ r,^ ( vx \rx y. 

9u\^ui 1 ^u) Su\^u 1 

Zu’" - Zf'^ 



) du — 

[ - 


t 

^x' 

) / ^x.i 
\^u 

-Zf)Y 




du 


T T 

+ I {gliX^, ir, - gliK^Yu, du-I z^’^ - z:'-* dWu. 

t t 

Thus, Girsanov’s theorem yields 


9\t,xn-9\t,xn 


< 


1 


eQ 

< 

where Q* is the probability measure given by 


xrX ryX 


^{\zy-z^'^\=o} d'^ I 


k,\X--X-\ + / ik^\x--Xl\ + k,\Y:i-Y:^\] du\Et 


dQ^ J gi{X-,Y-,Z-) - gliX-',Y-',Z-') 


dP 


= £ 


I yX,l yX',l 

\ ^ ^ 

By (A4’) and boundedness of Y^ and Y^' is well defined. Since by GronwalTs lemma we have 

T 

- xf I < (|Xf - xf\ + A :2 I \Y^ - Y^'\du)e^-'^, s G [t,r], 


^{|zS-*-zS'’Vo} ■ ^ 


it holds 

9\t,Xf)-9\t,xf) 


< 


+ Tk 3 )\X^ - Xf\ + Pk^Y k 2 k 5 e^^^) / \9{u, X^) - 9{u, X^')\du \ Pt 
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Hence, 


oHt,xn 


eKt,xf) 


< ut where ut is the solution of the ODE 


ut 


T 

+ Tk3)l\Xf - Xf\ + {k2k3Te^^^ + ki + j lus du 

t 


which is given by 


Ut = + Tk 3 )l\Xf - Xf'l exp + h + k 2 k 5 e'^^^)l{T - f)) . 


Thus, 


e{t,xn-0it,xn <K3 \x--x, 




with it's := + Tfcs)/ exp ((fc 2 A: 3 Te^i^ + ki + /c 2 A: 5 e^i^)/r) which show that 0 is a 

Lipschitz function and the Lipschitz coefficient does not depend on the bound M of Z. 

Let Ckxkd be the constant Ck^x,i,d with fcs replaced by K 3 and put N = [T/Ck^x^i^d\^ where 
[aj denotes the integer part of a, and ti := iCk,x,i,d, i = 0,... ,N and fAr+i = T. Since 
fi < Ck^x,i,d> by Theorem 2.1 the FBSDE 


iXt = x + /q bu{Xu,Yu) du + /g au dWu 

\Yt = gu{Xu,Yu,Zu)du- ZudWu, t E [ 0 ,fi] 

admits a unique solution {X^,Y^, Z^) such that \Z^\ < M with M = A\2K3\f^ for all 
t E [ 0 ,fi]. Therefore, {X^ ,Y^, Z^)l]^Q,ti\ — (X ,Y, Z)l[Q^t^y Similarly, we obtain a family 
y*, Z*) of solutions of the FBSDEs 


(Xt = Xu., + bu{Xu, y^) du + au dWu 
jy* = e{ti,Xu) + gu{Xu, Yu, Zu) du - Zu dWu, t E [ti-i,ti] 
such that {X\Y\ Z^)liu_,,u] = (X, Y, Z)l[u_,,u]> i = l,...,N + 1. Define 


N+l 

x-='E 


,h]> 


N+l 


N+l 

and Z:=Y. 


2 = 1 2 = 1 2 = 1 

Then, (X, Y, Z) E x x is fhe unique solution of the FBSDE (1.1) 

satisfying \Zt\ < M for alH E [0, T]. In fact, it is clear that (X, y, Z) E x x 

^cxD(]^Zxd) ^ finite sum of elements of the same space. Let t E [0, T] and f = 1,..., X + 1 
such that t E [fj-i, fj]. We have 


X + 





I tjAt tjAt \ 

J bu{Xi)du+ J GudWu 

Vi-1 ij-l / 


= XI = Xt 
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and 


T 


T 



t t 


That is, {X, Y, Z) satisfies Equation (1.1). 



□ 


4.3 Proof of Proposition 2.4 

By Theorem 2.2 the FBSDE 


Xt = X + bs{X, X)ds + /o* asdWs 

Yt = UXt) + It gs {Xs,Y„Zs)ds- ZsdWs , t G [0, T] 


has a unique global solution {X,Y,Z) G x x such that \Zt\ < M 

for some constant M > 0. Eet Xt = Yt = T^^X, Zt = T~^Z, we obtain that (X, X, Z) G 
X 5^(]R^) X such that \Zt\ < M for some constant M > 0. Moreover, 

(X, X, Z) satisfies fhe EBSDE (1.1). This complefes fhe proof. 


□ 
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